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Abstract. We consider equations of nonlinear Schrodinger type augmented 
by nonlinear damping terms. We show that nonlinear damping prevents finite 
time blow-up in several situations, which we describe. We also prove that the 
presence of a quadratic confinement in at least one spatial direction is sufficient 
to ensure that the solution of our model decays to zero for large time. In the 
case without external potential we prove that the solution may not go to zero 
for large time due to (non-trivial) scattering. 



1. Introduction 

Wc consider, for a > and A G R, the following class of damped nonlinear 
Schrodinger type equations (NLS): 

(1.1) ( + - + ^l-l"^- - ^ > 0' - ^ I^'^' 

where S denotes the energy space associated to the harmonic oscillator, i.e. 

E = {/ e H\^% X ^ \x\f{x) e L2(R'i)} , 
equipped with the following norm: 

ll^lls := \\u\\l^- + I|Vu||l2 + ||a;w||L2. 

In the following, we allow uq to be arbitrarily large within E, i.e. we shall not be 
concerned with solutions corresponding to small initial data. We make the following 
standard assumption on the nonlinearities: 

where (d — 2)+ denotes the positive part. Thus, if d ^ 2, we impose no size restric- 
tion on (Ti, (72 > 0. For d ^ 3 the above assumption ensures that both nonlinearities 
are If^-subcritical. The external potential V is supposed to be harmonic (or zero), 
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As we shall indicate below, all of our results can be generalized to the case of 
potentials V{x) > 0, growing at most quadratically at infinity. In the case a = 0, 
it is well known that (jl.ip is a Hamiltonian equation (its mass and energy are 



RC was supported by the French ANR projects SchEq (ANR-12-JS01-0005-01) and BECASIM 
(ANR-12-MONU-0007-04). CS acknowledges support by the NSF through grant no. DMS- 
1161580. 



2 



P. ANTONELLI, R. CARLES, AND C. SPARBER 



conserved); see e.g. [101 Hlj. In fact, the Hamiltonian counterpart of our model, 
i.e. with a £ iR (yielding a nonlinear Schrodinger equation with combined 

power-law nonlinearitics) has been studied in |20j . In the present case a > 0, the 
last term in (jl.ip is dissipative, which is the reason why we consider non-negative 
times only. Indeed, the dissipative nature of can easily be seen from the fact 
that the local conservation law for the particle density p = |up is augmented as 
follows: 

(1.2) dtp + divJ = ~2ap''^+\ 

where, as usual J ~ Im('uVu) denotes the current density. For a > the right hand 
side describes a nonlinear damping mechanism for the density. 

Equations of the form (jl.ip arise as phenomenological models in different areas 
of Physics. For example, in nonlinear optics, equation (|l.ip with V = models the 
propagation of a laser pulse within an optical fiber (d = 1) under the influence of 
additional multi-photon absorption processes, see, e.g., [5l|T2]. Another application 
arises from quantum mechanics, where NLS type models arise in the description 
of Bosc-Einstcin condensates in harmonic traps (which are experimentally required 
to produce these condensates). In this context the nonlinear damping is a model 
for the reduction of the condensate wave function through higher order particle- 
interactions, cf. [nil]. 

From a mathematically point of view, NLS type equation with nonlinear damp- 
ing terms have been studied in |13| as a possibility to continue the solution of NLS 
beyond the point of finite-time blow-up (see also |18j for an earlier study in this 
direction). In [3], global existence for the particular case ai = 1 and 1 < (72 ^ 2 
in dimensions d ^ 3 has been studied. Notice that in d = 3, this allows to take 
into account an if^-critical damping term (i.e. a quintic nonlinearity with 02 = 2). 
In [llj the particular case of a mass critical nonlinearity ai = 2/d and V = Q 
has been studied. In there, global in-time existence of solutions is established if 
(72 ^ 2/(i and it is claimed that finite time blow-up in the log- log regime occurs if 
(72 < 2/(i. A more complete understanding of the possibility of finite time blow-up 
remains an open problem, however (numerical simulations can be found in |18ill3) ). 

In the following, we shall develop a more systematic study of NLS type equations 
with nonlinear damping, generalizing the results mentioned above in several aspects: 

• We extend the results of global well-posedness to the case of general (energy- 
subcritical) nonlinearitics. 

• We prove that in the case without external potential, the solution is asymp- 
totically close to the solution of the free equation for t +00, i.e. we 
establish scattering for positive times. 

• We show that in the case where a quadratic external confinement is present 
in at least one spatial dimension, the norm of solution vanishes asymp- 
totically with a certain (not necessarily sharp) rate. (This result corrects a 
mistake in the proof of Corollary 2.5 in [5]). 

• We compare the results above to the ones which can be obtained for x G A/, 
a compact manifold without boundaries. 

Theorem 1.1. Let d ^ 1, a > 0, uji,. . . ,ujd ^ 0, and uq £ T,. Then, the Cauchy 
problem ()1.1|) has a unique solution u G C(R+;I]) in either one of the following 
cases: 
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(1) A ^ (dejocusing nonlinearity) and < cri,CT2 < 2/((i— 2)+; 

(2) A < (focusing nonlinearity) with, either 

(a) < (Ti < 2/d and < (72 < 2/(d - 2)+, or 

(b) 2/d ^ <Ti < (72 < 2/((i - 2)+, or 

(c) cTi = (T2 = 2/d, or 

(d) 2/(i < <Ti = (72 =: (T < 2/((i — 2)+ anii a ^ min [a, y/a) |A|. 

// in addition there exists at least one j such that ujj > 0, then u vanishes asymp- 
totically as t ^ +00; 

lim ||u(i)||i2(R.d) = 0. 

//, however, coj ~ /or al/ j, u may behave like the free evolution of a non-trivial 
asymptotic state u+ G L^(R'^); \\u{t) — e's^u+H^n-j^d) — > as t ^ +00. 

The paper is organized as follows: In the next section we establish the local (in- 
time) well-posedness of our model in the energy space. In Section [31 we extend this 
to global in-time wcU-poscdncss using a modified energy functional. The long time 
behavior and the possible extinction of solutions is studied in Section H) Finally, 
we briefly discuss the case of compact manifolds in the appendix. 

2. Basic properties of the Cauchy problem 

In this section we shall show that is locally wcU-posed for any uq G S and 
we also establish a blow-up alternative. 

2.1. Local well-posedness. We denote by U{t) = e~^^^ , the Schrodinger group 
generated hy H = ^1^^ + ^ ■ We first recall the standard Strichartz estimates (see 

e.g. m)- 

Definition 2.1. A pair {q,r) is admissible if 2 ^ r < (2 ^ r ^ 00 if d = 1, 

2^r<oo if d ~ 2) and 

Proposition 2.2 (Strichartz estimates). Let T > 0. There exists 77 > such that 
the following holds: 

(1) For any admissible pair {q,r), there exists Cq such that 

\\ui-MLH[a.v];Ln ^ c^Ml^, yip e l^cr''). 

(2) For s e R, denote 

Ds{F){t,x)= U{t-T)F{T,x)dT. 

For all admissible pairs (gi,ri) and (92, '"2); there exists C = Cq^^q^ independent of 
s S R such that 

(2.1) 

for all F € LP^I; L'^'^) andOs^S i^r]. 

(3) In the case without potential, V = 0, the above results remain true with 77 = 00. 
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Proposition 2.3 (Local existence). Let \,a E R, oji, . . . ,ujd ^0, and (Tj > with 
(7j < 2/{d — 2) if d ^ 3. For all uq G S, there exists T and a unique solution u of 
(|l.ip . such that 

u,Vu,xuEC{[Q,T]■L'{^R!^))C^ p| ([0, T]; L2-.+2(Rd)) . 

Proof. We present the main steps of tlie classical argument, which can be found for 
instance in [TOl [M] in the case V = Q (see also [8] in the presence of a potential) . 
Duhamel's formulation for reads 

(2.2) u{t) = U{t)uQ-i\ [ U{t-T){\uf'''u){T)dT-a [ C/(t-r) (jup'^^w) (r)dr. 



Denote the right hand side by $(it)(t). Proposition 12.31 follows from a fixed point 
argument in a ball of the space 



2crj+2 



= < u e C([0, T]; S) ; a;w, Vm € p| Lj,'"^ L 

[ J = l,2 

where L^L*^ stands for L9([0, T]; L''(R'*)). For j = 1,2, introduce the Lebesgue 
exponents 

/ , 4cr, +4 ^ 2cr,(2cr, +2) 

^^-^^^-^^ ^ '^^-^ ^ ^^- 2-(/2)4 - 

Then {qj,rj) is admissible, and 

1 _ 2crj 1 1 _ 2aj 1 

r'j rj rj ' q'. 0^ qj' 

Proposition 12.21 and Holder inequality yield, for j = 1,2 

m^)\\LyL':nL-L^^C\\uoh^-+C ^ IWuf^'u] 

£=1,2 



= 1.2 



where C is independent of T ^ 77. We note that if G( ^ 21 d for = 1,2, then 
Ot^qi, and 



(2.4) ll<f(-)ll,^.^..nL^L^ < C\\u4l^+CY. T2..(i/..-i/,.)|| 



||2^^ + l 



If (Tf > 2/(i for £ = 1 or 2, using Sobolev embedding, 

£=1,2 

We have 

Vi>(u)(0 = [/(t)VMo - «A / C/(t - t)V (litp^^w) (T)dT 

Jo 

U{t - t)V (IuP'^^m) (T)dT -i ( U{t-T) ($(w)(t)VF(t)) dr. 
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We estimate the second term of the right hand side as above, and, for all admissible 
pairs {q,r), the new term is estimated by 
ft 



C/(t-r)($(u)(r)Vy(T))rfr 



where we have written an estimate which is valid in the more general case where 
V is at most quadratic (d^V £ L°°(R'^) for \a\ ^ 2). Similarly, to estimate x^{u), 
a new term appears, which is controlled by 

CT||V$(u)|Uooi.. 

Choosing T sufficiently small, one can then prove that $ maps a suitable ball in 
Xt into itself. Contraction for the norm || • is proved similarly, and one 

concludes by remarking that Xt equipped with this norm is complete. □ 

Remark 2.4. The above result can be extended to the case where the first assump- 
tion is replaced with u £ L°°(R+; for some s > 0, up to changing the 
application of Holder inequality in the proof. 

Remark 2.5 (Energy-critical damping). When d ^ 3, the case cr2 = 2/{d—2) could 
be considered, like in [3] for the case d = 3. This requires a different presentation 
in the proofs, which is the reason why this case is not studied here. 

Remark 2.6 (More general potentials). As suggested in the course of the proof. 
Proposition 12.31 remains valid if we assume more generally that V{x) is smooth, 
and at most quadratic, i.e. d^V £ L°°(R'^) for all |a| ^ 2. 

2.2. Basic a priori estimates and blow-up alternative. In order to extend 
the obtained local in-time solution to arbitrary time intervals, we shall derive several 
a priori estimates. In a first step, we recall (|1.2p to infer the following. 

Lemma 2.7 (Mass dissipation). The local in time solution u(t) e S satisfies 
d 

di' 

As a consequence, we have that u G L°°([0,T] x L^^R'')) n L2'^^+2([0, T] x R''). 

Proof. We multiply (jl.ip by u and integrate with respect to a; G R"*. Taking the 
real part, yields (|2.5p . which consequently implies 

d_ 

di' 

and thus ||it(t)||/^2 ^ ||-uo||l2, Vi G [0,r]. In addition, we can integrate (|2.5p with 
respect to t to infer 



(2.5) ^Ju{t)\\l, + 2a\\u{t)\\lZt^, =0, Vt G [0,r]. 



;\\u{t)\\l, ^ -2a\\u{t)\\lZt'^ 



2a / \\uit)\\%ltl. dt = lluolli. - \HT)\\l, ^ ||«olli2, 
Jo 

and thus u G L2<^^+2qo,T] x R'^). □ 

Remark 2.8 (Non-existence of steady states). An immediate consequence of (|2.5p 
is the non-existence of non-trivial steady states. In the Hamiltonian case (a = 0), 
they are found by inserting the ansatz u(t, x) = '0(a;)e'''* with /i G R into and 
study the resulting elliptic equation for ■)/). In our case, (|2.5|) together with the fact 
that for stationary states |u(t,x)p = |'0(a;)p, immediately implies that ijj ~Q. 
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Proposition 12.31 and Lemma [2?7l allow us to infer the following blow-up alterna- 
tive. 

Corollary 2.9 (Blow-up alternative). Let A e R, a,aji, . . . ,0;^ ^ 0, ctj > with 
<7j < 2/{d — 2) if d ^ 3, and uq G S. Either the solution to (jl.ip exists for all t ^ 0, 
i.e. 



(2.6) u, Vk, xueC (R+; L2(R<i)) n f| L.^'J^ (R+; L'-^+^{R'')) , 
or there exists T > 0, such that 

< 

In the case (Ti = (72 = f ; or m the fully mass-subcritical case ai,a-2 < 2/d, the 
solution is global, that is, (|2.6p is satisfied. 

Proof. Let M > 0. Lemma [2771 shows that 1 1— >■ ||w(i)|||2 is non-increasing function. 
Thus, the only obstruction to well-posedness on [0, M] is the existence of a time 
< T < i\/ such that 

\\xu{t)\\L2 + \\Wu{t)\\L2 ^+00. 

< 

As long as w G C{[0,t]; S), we have 



dt 



— / xAu{t,x)\ dx = 2Ke / Xju{t,x)dtu{t,x)dx = 21m / x.ju{t,x)idtu{t,x)dx 



Im 



x'^ju{t,x)Au{t,x) - a J x]\u{t,x)\^"^+'^dx 



2Im / Xju{t,x)dju(t,x) ~ a i a;||M(t, a;)p°'^+^da; 



^ 2\\xu{t)\\LA\^u{t)\\ 



where we have used Cauchy-Schwarz inequality and the assumption a ^ 0. Suppose 
u € L°°([0, T]; iJ^). Then the above estimate and Gronwall's lemma show that 
xu G L°°([0, T]; i^), hence a contradiction. Hence the first part of the corollary. 
The second part follows from the standard criterion for L^-critical problems (see 
e.g. [10]): if the solution does not satisfy <\2.Q\ . then there exists T > such that 

l|w(t)!|^t+vdrft = oo. 

This is also a direct consequence of the proof of Proposition 12.31 (if a\ = ai = 2ld, 
then Qj = qj). Lemma [2.71 rules out this possibility; this point has already been 
noticed in [TT| . 

Finally, if cri , (72 < 2/d, the standard argument given initially in [21] , following 
again from Proposition [231 shows that (|2.6p follows if m G L°°(R+; L^(R'')), which 
in turn is a direct consequence of Lemma 12.71 □ 
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3. Global well-posedness 



The a-priori bounds obtained in Seetion 12.21 are not sufficient to infer giobal 
weli-posedness for crj ^ 2/d (unless ui = <T2 = 2/d). In order to obtain further 
a-priori estimates, one possible approach would be to follow jllj . where the author 
studies the time-evolution of ||Vit(i)||L2 and shows that for CT2 > cti = 2/d one can 
obtain a bound of the form 

\\S/u{t)\\L2 s; \\Vu„\\l2 

for some C > depending on the involved parameters a, A, (T2- Indeed, an analogous 
result could also be obtained in our, more general, situation. However, we shall 
rather follow the approach of [3] based on a modified energy functional E{t) which 
will allow us to infer (under certain conditions) uniform in-time bounds on different 
quantities involving u{t). 

3.1. Bounds on a modified energy functional. In the following, we denote 
E{t) = \\\Vu{t)\\l, + f V{x)\u{t,x)\^dx 

(3.1) X r r 



\u{t,x)r'+'dx + K / \u{t,x)r^+Ux, 



(Tl+1 

for some k > (to be made precise below). Clearly, E is well defined on [0,r] for 
any u G C([0, T]; S), by Sobolev embedding. Even though this energy functional is 
not conserved, we shall prove that it is uniformly bounded in time, provided that 
some assumptions on the involved parameters hold true. 

Proposition 3.1 (Energy bound). Let Q < k < and assume that 

(1) Either A ^ 0, 

(2) Or X < and a2 > (7i. 

Then there exists a C ~ C{\\uq\\]^2) ^ such that: 

E{t)^E{0)+Ci\\uo\\L2), Vt€[0,r], 
where T > is the existence time obtained in Proposition \2.3\ 



Proof. We first assume that u{t) is sufficiently regular an decaying so that all of 
the following formal manipulations can be carried out. Once the final result is 
established, a standard density argument allows to conclude that it also holds for 
weC([0,T];E). 

We compute the time-derivative of the energy functional (|3.ip . using equation 
(|1.1[) . which yields: 

(3.2) ^E{t) = a [ Re(uA?l) - K(cr2 + 1) / \u\'^''nm{uAu) dx 

-2a f Vix)\u\^^'+^dx-2aX f \u\^"'+^^'+^ dx 
-2aK(cr2 + l) / \u\'^^^+^dx. 

jR-i 

Consider the first term on the right hand side; since A|wp — 2Re('i2Au) -I- 2|Vwp, 
notice it can be rewritten, using integration by parts, 

I Iwp'^^ Re(uAu)da; = - /" iVw^da; - 2cr2 / \u\^''^\W\u\\^ dx. 

JR'' Jr-' Jr-' 
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Now, we rewrite the second term on the right hand side: 

/ lup""^ Im(itAu)dx = / lup'^^iv (lm(ftVu))da; = - / V\u\'^''^ ■lm{uVu)dx. 

Here we use the polar factorisation introduced in [5] (see also [3]), to show the 
above integral equals 

(3.3) - 2(72 / lup"' Re(0Vu) • lm{$\7u)dx, 

where 4> is the polar factor related to u, 

|o iiu{t,x) = 0. 

This indeed can first be proved by replacing (f) with 

y/\u{t,x)\^ 

and then passing to the limit e — > in H^, as in [5], [5]. Let us rewrite p.3p as 
0-2 y" liip'^'l Re(0Vw) -Im(0Vu)f - lu]^"^ ^|Rc((^Vu)p + | Im((^Vu)p^da; 

= (72/ |w|2'"^|Re(^Vu)-Im(0Vw)|^dx-cr2 / luf^^l'^^ufdx. 

The last equality follows from the identity 

|Vwp = I Rc((^Vu)p + I Im(0Vu)p, a.e. in R'', 

see formulas (30) in [5] and (5.15) in [3]. Hence, by resuming the second term on 
the right hand side of p.2p is equal to 



-K(cr^ + (72) / \uf''^\Rc{(j)\7u)~lm{(j)\7u)fdx + K{a2+(J2) [ luf^'lVu^dx. 



In summary this yields: 
d 



E{t) = - (a - K{al + (72)) / litp'"^ \Vu\^dx - 2a(72 / |V|u 



|2 



-k((7| + (72) / |u|^'^^|Re((/)Vu)-Im((/iVu)|^da;-2a / V"(2:)|u|^'^^+^ da; 

-2aX f |u|2'^=+2<Ti+2^^_2aK((72 + l) / \u\'^'''+'^ dx. 

Under the assumption < k < ^3^^, and if A ^ (dcfocusing case), all the terms 
on the right hand side are non-positive, and we infer that i? is a non-increasing 
function: E{t) < E{0), yt e [0,T]. 

If A < (focusing case), however, the term involving the 7^20-1 +2cr2 +2 ^-^qj-^-^ jg 
positive. But since (72 > (7i by assumption, wc can interpolate this term using: 

||u||l2<ti+2<,2+2 < ||'"|Il4<t2 + 2 ||u||^24 + 2 J 

with 

1 e , . ^_ (7l(2(72 + l) 



0'l+(72 + l 2ct2 + 1 (72 + 1' (72(cti -I- (72 + 1) 



e (0,1). 
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Denoting 7 = cri/fT2, this implies that 

II ||2ai+2<T2+2 ^ II ||7{4<T2+2) II ||{1-7)(2t2+2) < p||,. ||4<T2+2 , I ||,/||2'^2+2 

II "II ^2<Ti + 2,T2+2 ^ II "11^4,^2 + 2 II "II ^2^2+2 ^ ^ II "II ^4,72 + 2 ^ ^7/(1-7) L2<t2+2 I 

where we have used Young inequaUty. In summary this yields 

^^E{t) ^ - 2a(.(a2 + 1) - |A|.) J^^ |u|^-+^ dx + J^^ dx. 

For < e <C 1, the coefficient in front of the first term is negative, which allows to 
conclude, after an integration with respect to time, that: 

yt e [0,T] : E{t) E{0) + C [ [ \u\^'''+'^ dx dt < E{0) + C{\\uo\\l^), 

since we have u € L'^'^^+'^{[0, T] x R'') from Lemma [2Jl □ 

We are now in the position to prove global wcll-posedness of (jl.ll) under various 
conditions on the parameters. For the sake of a simpler presentation, we shall treat 
the case cti = (72 separately, see Section 13.31 below. 

3.2. The case cri ^ (J2- In the following we shall consider two different powers and 
impose the following assumption. 

Assumption 3.2 (Nonlinearity). Let a\ ^ <J2 and, in addition: 

(1) Defocusing case. //A ^ 0, we assume 

2 

< (71,(72 < ^^—^ (0 < (71,(72 i/d < 2). 

(2) Focusing case. If X < 0, we assume 

2 ■? 
• Either, < (7i < — and < (72 < 



d ^ {d-2)+' 

Remark 3.3. Assumption 13.21 can be understood as follows: If without damping 
(a = 0), the solution of is global in time, then the strength of the nonlinear 
damping plays no role. On the other hand, if finite time blow-up may occur in the 
Hamiltonian case (i.e., A < and cti ^ 2/(i), then the damping is assumed to be 
stronger than the attractive interaction, in terms of the power (72 . 

Theorem 3.4 (Global existence I). Let d ^ 1, a > 0, and oji,. . . ,uJd ^ 0. Under 
Assumption lS.Sl for any uq G S, (jl.ip has a unique solution 

4<Tj+4 

C(R+;E)nL°°(R+;i7i)nL2<^^+2(R+xR'^)n f] {R+; L^''^+^ {R'^)) . 

J = l,2 

Moreover, if iOj > for all j S {1, . . . , li}, then we also have u e L°"(R+; E). 

Proof. In the defocusing situation A > 0, Proposition 13.11 immediately implies a 
uniform (in-time) bound on ||VM(t)||^2 ^ E!(0), since E(t) is the sum of four non- 
negative terms. In view of the blow-up alternative, stated in Corollarv l2.9[ we thus 
infer global well-posedness. 

For the focusing situation A < 0, we note that the mass-subcritical case (7i, (72 < 
2/d, has already been dealt with in Corollarv 12.91 (including ai = (72 = 2/d). The 
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moment ai ^ 2/d, we require a2 > o\ which aUows us to interpolate the potential 
energy of the attractive nonlinearity similarly to the calculation above, namely: 

for any e > and an appropriately chosen < /? < 1. Now fix e = 2k{g\ + 1)/|A| 
to obtain 

for some C — C{X,K,(7j) > 0. In view of CoroUarv 12.91 and Proposition 13.11 this 
yields ||Vm(<)||^2 ^ E{0) + llwollia and we are done. □ 

3.3. The case cti = o'2- In the case where cti = (72 = cr, (jl.ip simplifies to 

(3.4) idtu+^Au = V{x)u + {X-ia)\u\'^''u, U|i^o = e S. 

Formally, this model can be considered as the diffusionless limit of the (generalized) 
complex Ginzburg-Landau equation. Concerning global existence of (|3.4I) , we shall 
only be interested in the L^-supercritical situation (cr > 2/c?), in view of the last 
assertion in Corollary [ 



Theorem 3.5 (Global existence II). Let d ^ I, a > 0, a > 2/d, anduj\, . . . ,uJd ^ 0. 
Assume that 

(1) Either A ^ 0, 

(2) Or A < anrf a ^ min (cr, ^) |A|. 
Then, for all uq G S, (|3.4p has a unique solution 

u e C(R+;I])nL°°(R+;i/i)nL2<^+2(R+ x R'^) (R+; i2'^+2(R'^)) . 

Moreover, if ujj > for all j G {1, . . . , d}, then we also have u G i°°(R+; E). 

Proof. The defocusing case A 0, can be treated as in the proof of Theorem 13.41 
above. It therefore remains to consider the case A < and cr > 2/d. Note that we 
cannot invoke Proposition 13.11 which requires cri < cr2. Instead (following an idea 
in [3]), we consider the linear energy functional: 

EUt) = ]-\\Vu{t)\\l,+ f Vix)\u{t,x)\-'dx. 

Obviously, it suffices to prove that Eiin{t) is uniformly bounded in time. Differen- 
tiating E\in{t) and using Equation p.4p yields 



(3.5) 



Elicit) = a / lu]^" Re{uAu) dx + |A| / \u\^" Im(wAu) dx 
dt J^d J^d 



-2a / V{x)\u\^^+^dx. 

Using the same arguments as in the proof of Proposition 13. 1[ we infer 

l EUt)^^ia^\XW) f \u\^"\Vu\^dx^2aa f \u\^- Mu\\^ 
at j^d jjid 

-\\\a I |up'^|Re(0Vu)-Im(^Vu)pda;-2a / V[x)\u\^''-^'^ dx. 

Jr^ jRd 
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Thus, a a\X\, we obtain E\in{t) ^ i?iin(0) < oo, for all t e [0,T]. On the other 
hand, from p.Sp we have 



dt Jjf^d Jj^d 

+ 2a\X\ [ \u\^'^\7\u\-lm{4>Vu)dx-2a [ V {x)\u\^'^+^ dx . 

By using Caucliy-Schwarz and then Young inequality, we see that the third term 
in the right hand side is bounded by 



hence if IAI-^ct < a, then again we have -Eiin(i) ^ i5iin(0) < oo, for all t e [0,T]. 
This establishes global well-posedness of ()3.4p . in the second case. □ 



Remark 3.6. The second case of Theorem 13.51 shows that if the damping is of the 
same strength (in terms of its power) as an attractive interaction nonlinearity, the 
relative size of the respective coefficients starts to play a role. Note that the larger 
the dimension d ^ 1, the smaller a can be chosen to ensure global existence. We 
remark that the numerical experiments presented in |13j always consider a ^ |A| 
and thus they do not show the aforementioned possibility for global existence. 

At this point, we have proved global existence in all the cases listed in Theo- 
rem 11.11 Wc now turn to the large time behavior, in cases where the solution is 
defined for all t ^ 0. 



4. Large time behavior 



In view of the dissipation equation (jl.2p . the damping is expected to have a 
significant influence on the long time behavior of solutions to (jl.ip . Indeed, if we 
consider, for a moment, the case of x-indcpendcnt solutions, then ()1.2|) simplifies 
to the following ordinary differential equation 

dtp = -2ap''-^+\ p|t^o = Po:=|wo|', 
the solution of which is given by 

Po 



(4.1) Pit) = 



(l + 2atp;^^)i/-2 



Thus, one might expect the solution to vanish like ||M(t)|||^2 = 0{t~^^'^^), as i — ?> 
+00. We shall see below, however, that this rather naive argument does not yield 
the correct long time behavior of u in the case of where ^ = 0. The idea is that 
the dispersion due to the Laplacian may prevent the damping from taking the wave 
function u to zero. On the other hand, in the presence of a confining potential (i.e., 
ujj > for some j), there is no (large time) dispersion, and the ODE mechanism 
remains, at least qualitatively. 

4.1. Asymptotic extinction with full confinement. In this subsection, we con- 
sider the case of a fully confining potential, in the sense that we suppose ojj > 
for all j. We start with a simple estimate, which can be viewed as a dual version 
of Nash inequality. 
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Lemma 4.1 (Localization). Let 1. For all p ^ 2, there exists C such that for 
all f e S{R'^), 

Proof. For R > 0, write 

ll/llL2(Rci) ^ \\f\\L^{\x\<R) + \\f\\L^{\x\>R) 

where we have apphed Holder's inequality for the first term, and simply multiplied 
and divided by |a:;| for the second term. Both terms on the right hand side have the 
same order of magnitude if R^+'^(^/^~^/p'> = ||a::/||L2/||/j|ip. Using this value of R 
yields the result. □ 

Recall that in the case with confining potential {uij > for all j G {1, . . . , d}), 
the solution satisfies u G L°°(R+; E). This can be used to obtain an estimate for 
the time-decay of the solution. 

Proposition 4.2 (Asymptotic extinction I). Let d ^ 1, a > 0, uji,...,u!d > 0, 
and uq GT,. In either of the cases mentioned in Theorem \l.l[ the solution to (jl.ip 
satisfies u G L°°(R+; S) and there exists C > such that 

||u(t)||^2 Cr'w^, Vi ^ 1. 



Proof. In view of Lemma [4.11 with p = 2(72 + 2, we have, since u e i°°(R_|-; S) by 
assumption. 

Along with (|2.5p . this yields 



-\\uit)\\i.+2aC\\uit)\\^./ ^0. 
Therefore, y{t) = ||u(t)|||2 satisfies a differential inequality of the form 

m + cyitr ^0, p^i + ^^a^. 

The comparison with the ordinary differential equation yields y{t) = 0(i~^/(P~^)) 
for i ^ 1, since p > 1, hence the result. □ 



We shall see in Section 14.21 that the presence of at least some confinement is 
essential for the time-decay of the norm. 

4.2. Non-vanishing solutions, hi the following, we assume that there is no con- 
fining potential, i.e. = in p.ip . Then the solution must not be expected to 
vanish as i — > +oo, as shown by the following result. 

Proposition 4.3 (Scattering without potential). Let d ^ 1, A,o G R and 2/d ^ 

fi,cr2 ^ 2/((i — 2). Let R > 0. There exists T depending only on d, \\\, \a\ 
and R such that for all u+ G H^^Rf') with \\u+\\hi ^ R, there exists a solution 
u e C([r,oo);Hi) to 

idtU + ^A?! = A|up'"iu - ia\u\'^'"^u, 
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such that 

lim \\u{t) - e'^^u+\\Hi(Rd) = 0. 

Sketch of the proof. This result follows from [TS] (see also [TH Proposition 3.1] for 
a simplification), and stems from a fixed point argument applied to Duhamcl's 
formulation of the above problem: 

/•oo />oo 

u{t) = U{t)u+ + iX U{t - s) (H^^'u) {s)ds + a U{t-s) {\u\^''^u) {s)ds. 



Denote by the right hand side, and Ui^cc{t) = e'2^u+. Let rj = 2a j + 2 

and qj > 2 be such that {qj,2aj + 2) is admissible, j = 1, 2, like in (|2.3p . With the 
notation L^Y = i'^([r, oo); F), we introduce: 

=^ 2j|it+j|Hi , ||u||^9,-^2^j+2 ^ 2 ||itf„o||^;^j^2„,+2 , j = l,2|, 

where Cq. is given by Proposition 12.21 Resume the numerology of the proof of 
Proposition [231 we have 

1 _ 1 2(Tj 11 2(7 

rj rj q, 9^ 

where qj ^ 6j < oo since 2/d ^ aj < 2/{d— 2). For u £ Xt, Strichartz estimates 
and Holder inequality yield: 

||$(u)|l4.,^i,2.,+2 Cq^\\u+\\H^+c J2 (lll"l'"'"IL<^-;, + lll"l'"'^"IL'Ji.O 

£=1,2 ^ 
^=1,2 

for ?7£ = (/f/^f £ (0, 1]. Sobolcv embedding and the definition of Xt then imply: 

1=1,2 

We have similarly 

£=1,2 

1=1,2 

From Strichartz estimates, ufrco G L'^^ (R;£'^), so 

||ufroc||^;y ^2^j+2 ^0 as r ^ +00. 
Since r]£ > 0, we infer that $ sends Xt to itself, for T sufficiently large. 
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We have also, for U2,wi e Xt'- 



1=1,2 



< \\uh^cfjZ''l.J\u+fH!'-^ '"'^ ll"2 -uiI1l:«l., ■ 



£=1,2 



Up to choosing T larger, $ is a contraction on Xt, equipped with the L^L"^^ n 
L^^L''2-norni. Since this space is complete (see e.g. [TOl Section 4.4]), the proposi- 
tion follows from the standard fixed point argument. □ 



Proposition 14.31 rules out the asymptotic extinction of the solution, since, by 
invoking the triangle inequality, we directly infer 



lim 



|i2 ^ lim ||e'3^u+|U2 - ||M(t)-e*5^w+||i2 



: lim ||e'2^w+||2^2 

t-^-\-QO 



due to mass conservation. The solution of the damped equation therefore does not 
decay to zero, due to the possibility of radiation escaping to infinity. The latter is 
no longer true if we consider (jl.ip on a compact manifold, instead of R'^ (see the 
appendix), or in the presence of a partially confining potential on R'', as we show 
in the next subsection. 

4.3. Asymptotic extinction with partial confinement. Suppose d ^ 2, and 
decompose 



X = {x',x") e R'' X R' 



d-k 



with > if 1 ^ j ^ k and luj = if j ^ k + I. Consider the linear equation 



(4.2) 



1 A 



In this case, a generalized Mehler formula is available (see, e.g., [17] )j yielding the 
linear solution for t > 0, in the following form 



1 1 /" ^.^ ^ y^.^—i ( — ^COSt — Xjyjj - -y' 

uiin{t,x) = — ■ — -jTF / e"" V J e' 5* uo{y)dy, 

(2z7rsini)2 [2nTt) 2 jRd 

where x = (cci, x') and the square roots are defined in such a way that Maslov index 
is taken into account. In the case where uq is a tensor product, uo(a;) = f{x')g{x"), 
we have 

uUt,x) = ('e~K-^''-'' + ^)/(x')') ® (e'^^--'=5(x")) ■ 



Even though the evolution of / leads to no global dispersion (the solution is periodic 
in time), the standard dispersive properties of the free Schrodingcr group allow to 
write that, for 2 < p < 00, it holds: 



klin(t)|l 



LP(R'') 

O (t-('^-fc)(i/2-i/p)) as t ^ +00. 



LP(R''-'=) 
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For a general uo G S, it is possible to recover such a decay rate by invoking suitable 
vector fields. In view of the relations (see e.g. [8]) 

A{t) = e-"^ (zV,.,x") e''"f = (- (sini) x' f + i (cost) V,,/, (.t" + itV ^n) f) 

i (cOsOe-'^'^'^'V^- J^g»l^tanty^ ^^^^.i^li^^^^ 

= e-**^xe'*^/ = ((cost) x' f + i (sint) V,-/, (x" + iiV,») /) 

i (sint)e'^™**V^- j^e"*^™'*/^ ,ite'^V^" J^e"'^/ 

the first equality shows that A and B commute with the linear group e~'*^, so 
||^(t)'"iin||L2 and ||_B(t)iiiin||L2 arc conserved quantities, and the last equality, com- 
bined with Gagliardo-Nircnberg inequality, yields 

C 2 

hlin(t)||LP(R<i) < ^(rf-fc)(l/2-l/p) H"0||s, 2^p< J^—^- 

In view of the previous subsection, the effect of nonlinear damping terms, in combi- 
nation with partial confinement, does not seem obvious, because without damping, 
there may be scattering for the nonlinear analogue of (|4.2p (provided that the 
nonlinearity in short range in R''"'''). 

The present situation can also be compared to the framework of [IGj and [22] . 
There, cubic or quintic Schrodinger equation is considered on a product space of the 
form R'^~'^ X M'^, where M*^ is a fc-dimensional compact manifold, and scattering 
is established at least for small data. Let us describe more precisely the results in 
[16j : The quintic defocusing equation is considered on R x T^, a situation which 
may be compared to the case of a partial harmonic confinement in the last two 
spatial directions: 

i9tu -f Arxt^u = ; u\t=o = uo- 

Z. Hani and B. Pausader prove that even though global dispersion is not possible 
in the last two components, dispersion associated to the first component is enough 
to ensure scattering: assuming that ||wo||j?i(RxT2) is sufficiently small (a condition 
which is conjectured to be unnecessary in [l^), then there exists v+ G i/^(R x T^) 
such that 

\\u{t) ^ e'*'^^>^-r-'v+\\ Hi (RxT^) > 0. 

In [22], a similar result is proved for cubic (focusing or defocusing) Schrodinger 
equation on R" x M^, n ^ 2, in a more cumbersome function space. In view of 
Proposition 14.31 one might believe that adding a nonlinear damping is not enough 
to drive the solution to zero (in i^) for large time. 

We shall see that nonlinear damping always leads to asymptotic extinction in the 
case of a partial confinement, provided that the solution to is defined for all 
t ^ 0. The key to prove this assertion is an anisotropic counterpart of Lemma |4T] 

Lemma 4.4 (Anisotropic localization). Let d ^ 2 and 1 ^ k ^ d — 1. For all 

(72 > 0, there exists C such that for all f G iS(R''), 

llfl|2 r^ll f||(2<^2+2)/((fc+l)a2 + l)|| / .||(2fcfT2)/((fc+l)fT2-Hl) 

where x = {x\x") G R'' x R'^"''. 
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Proof. Let 

p(x')= / \f{x',x")\'dx" 

and write 



p{x')dx' = / p{x')dx' + / p{x')dx'. 

R'' J\x'\<R J\x'\>R 



The second integral can be estimated exactly like in the proof of Lemma 14.11 i.e. 

p{x')dx' / \x'\^p{x')dx' = R-^x'f\\l,^^,y 

x'\>R J\x'\>R. 

For the first integral we estimate (this step is analogous to Lemma 14. ip : 

\ l/('T2+l) 

p{x')dx' ^CkR'"''/^"'+'U I \pix')p+^dx'] 



\x'\<R. \J\x'\<R. 

In view of Minkowski's inequality, we have 



\pix')r+'dx'] ^ ( dx'. 

I J\x'\<R \JR'i-'= / 



|2-'|<-R 

Using (the same) Holder inequality, we infer 

1/(ct2 + 1) 

^ ^. DfcCT2/(cr2+l)ll /II 2 

2(R<i) 



(L . \fix',x")\'^'+^dx"y^''^^'\x' ^ C,i?^-^/(-^+l)||/||i2.2 + .(R. 



x'\<R yJii''-'' 
and for all i? > 0, 

ll/lli=(R.) < c',i?2'='^^/(-^+i)||/||i2.2+.(R.) + i?-'lk7lli.(R.). 

Optimizing in R yields the result. □ 

In view of Section [3] and of the proof of Proposition 14.21 we readily infer: 

Corollary 4.5 (Asymptotic extinction II). Let d^2, l^k^d~l, a>Q, 
LJi, . . . jLOk > 0, UJj — for j ^ k + 1, and uq G S. In either of the cases 
mentioned in Theorem the solution to satisfies u G L°°(R+; iJ^(R'')), 

x'u e L°°(R+; L^(R'')), and there exists C > such that 

||M(<)||i2 s$ cr T^^+fe, vt>i. 

Remark 4.6 (Comparison with the fully confining case). The decay rate is always 
weaker than in ProDOsition l4.2l Technically, this is due to the fact that in the proof 
of Lemma l4.4| we have used Minksowki's inequality and an extra Holder inequality, 
so we have lost more information than in the proof of Lemma 14.11 Note however 
that none of the rates in Proposition 14 . 2 1 or Corollarv l4.5l is claimed to be sharp. 

We will see in the appendix that in the geometry considered in [16] and [22] , the 
introduction of a nonlinear damping term also leads to asymptotic vanishing of the 
solution, in the same spirit as the above result. 
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Appendix A. Damped NLS on a compact manifold 

In the following wc shall consider x £ M, a d-dimcnsional compact Ricmannian 
manifold without boundary. A particular example is M = T'' = {'R/2ttZY, the 
d-dimensional torus. Since M is compact by assumption, we do not gain anything 
from the inclusion of a confining potential V. We thus set uij = 0, for all j = 1, . . . ,d 
and consider the equation 

^^^^ (idtu+^Au = X\u\'^'^^u-ia\u\'^'''-u, t^O,xeM, 

where A denotes the Laplace-Beltrami operator on M. 

The Hamiltonian analogue of (|A.2p has been studied in [7]. Obviously, we can- 
not expect the dispersive nature of the free Schrodinger group U{t) to hold on 
a compact manifold (a simple counterexample is given by the eigenfunctions of 
the Laplace-Beltrami operator on M). It turns out that this is true even locally 
in-time, see Remark 2.6 in [7|. The possibility of obtaining Strichartz type esti- 
mates on M is therefore severely restricted and any proof of a possible bound on 
l|f^(i)/l|L9(/;L'-(A/)) requires totally different techniques from those needed in the 
case M = R''. In view of this, we shall restrict ourselves to the situation d ^ 3, 
only, in which case the following local well-posedness result can be directly deduced 
from [7] (the case d = 1 is not treated in [7], as it is straightforward using the 
Sobolev embedding H^{M) ^ L°°{M), that is, without Strichartz estimates): 

Lemma A.l (Local well-posedness on compact manifolds). Let A, a e R, and 

assume that it holds 

(1) cri,(T2 > 0; if d ^ 2, and 

(2) < 0-1,(72 s: 1, i/d = 3. 

Then, for any uq G H^{M), there exists a time T > depending only on ||Mo||i/i(M) 
and a unique solution u G C([0, T]; iJ^(A/)) of (jA.2|) . 

Remark A. 2. In the particular case M = T'' several additional results are available. 
For example, in d = 1 local well-posedness in holds for nonlinearities which are 
smaller that quintic [B]. In higher dimensions, though, the situation seems to be 
more involved (at least if one seeks for strong solutions in the energy space). 

The proofs given in the case of R'' readily yield the following result: 

Lemma A. 3 (A-priori estimates). Let d ^ 3, A, a € R, cti, a2 > 0, with cti, CT2 ^ 1 
ifd ~ 3, and uq <E H^{M). On any time intervals I 3 such that u € C(/; H^(Al)), 
we have: 

(1) Mass dissipation: 

|||u(t)||i.+2a|l«(i)|li?.^-?.=0, yteL. 

(2) Control of the modified energy: for < k < , let 

E{t)^\\\Vu{t)\\l. + ^— [ \u{t,x)\'^^+^dx + K f \u{t,x)\'^^+^dx. 
2 CTi + 1 Jm Jm 

If either A ^ or A < and a-i > g\, there exists a C = C(||uo|1l2) 5^ 
such that: 

£;(t) «;-B(o) + c(||«o||l^), Vie/. 
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Proposition A. 4 (Decay rate on compact manifolds). Let d ^ 3, X,a E R, 
CTi,cr2 > 0, with (Ti, (72 ^ 1 */ = 3, and uo G H^{M). If either A or 
A < and 02 > <Ji, (IA.2|) possesses a unique global solution (in the future), 
u e C(R-)-; -ff^(-M)). In addition, we have 

\Ht)\\l2^ , Vt>0, 

and thus u vanishes asymptotically as t — > +00. 

Proof. Global existence for positive time stems directly from Lemma lA.ll and 
Lemma I A. 3 1 Since M is compact, Holder inequality yields 

We infer form the equation of mass dissipation, to get 
d 

di' 

The ODE mechanism sketched in Section 2] yields 



^Juml.+2a\Mr\\u{t)\\l7+'^0. 



\Ht)\\L2 ^ 



(l + 2a|Mh^|l^.o|li?)'^"^' 
hence the result. □ 

The rate established above is the one indicated by the naive ODE argument (|4.1I) . 
Obviously, the decay on M is faster than in the case of partial, or full confinement, 
but as underlined above, none of these rates is claimed to be sharp. 

Remark A. 5. The case of compact manifolds corresponds to the one found in nu- 
merical simulations (for example, M = T** in the case of spectral methods) , so the 
rate of Proposition IA.4I should be (at least) the one observed numerically. 

To conclude this paper, we examine the case shortly presented in Section 14. 3[ 
namely the case where the spatial domain is of the form M'^ x R''"'^, with k ^ 1. 

Proposition A. 6. Let k,n ^ 1, and M a compact manifold of dimension k. Let 
a > 0, A G R, (Ti, (72 > wq G S, and suppose that the following Cauchy problem 

^^2) H9tM+ ^Aa/xR"M = A|Mp'^iM-m|up'^^u, t ^ 0, x e M x R", 

has a solution u G C(R-|-;L^), satisfying 
d 
It 

Then u vanishes asymptotically as t +00 



u{t)\\l2 + 2a\\u{t)\\l%lt^, =0, yt^ 0. 



iWlli^ TT > Vi>0. 

(2at)'/"^|M|2' 



Proof. It suffices to notice that Lemma |4~41 is readily adapted to the present frame- 
work. Writing x = {x',x") £ M x R", and let 



p{x')^ I \f{x',x")\'dx". 

R" 
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Holder inequality yields 

p / f \ l/(<^2 + l) 

J p{x')dx'^\M\'''^^^'+^Uj \p{x')p+^dx'] 

In view of Minkowski's inequality, wc have 

/ f \ i/(t2+i) r / r \ i/(<^2+i) 

\pix')r+'dx'j ^ J^^ i^J^^ \fix',x")\'^^+'dx"j dx'. 

Using (the same) Holder inequality, we infer 

M \JR" / 

hence 

ll/lli^(MxR") ^ |Mp-/(- + l)|l/|ii2„.+2(MxR")- 

We can then resume the proof of Proposition I A. 41 with \M\ replaced by |Afp. □ 
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